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Abstract 

In this paper we present the effects produced by the temperature 
^ . in the renormalized vacuum expectation value of the zero-zero com- 

ponent of the energy-momentum tensor associated with massless left- 
handed spinor field in the pointlike global monopole spacetime. In or- 
der to develop this calculation we had to obtain the Euclidean thermal 
Green function in this background. Because the expression obtained 
for the thermal energy density cannot be expressed in a closed form, 
its explicit dependence on the temperature is not completely evident. 
So, in order to obtain concrete information about its thermal behavior, 
we develop a numerical analysis of our result in the high-temperature 
limit for specific values of the parameter a which codify the presence 
of the monopole. 
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1. Introduction 



The Euclidean Green function associated with a massless left-handed spin- 
1/2 field in the pointlike global monopole spacetime, Sf(x, x'), was presented 
recently in Ref. |J. There, Sf(x,x') was obtained applying the Dirac op- 
erator TpL, on the two-components bispinor Green function Gf{x,x'), which 
obeys a non-homogeneous second order differential equation similar to the 
scalar case. Using this procedure we can write 

S F (x,x')=f L G F (x,x'). (1) 

The massless left-handed spinor thermal Green function, S/3f{x,x'), can 
also be obtained from the respective thermal bispinor G/3f{x,x') in similar 
way. This fermionic Green function should be antiperiodic in the "Euclidean" 
time t with period (3 = 1/kbT, Kb being the Boltzman constant and T the 
absolute temperature. 

The gravitational effect produced by a global monopole can be approxi- 
mately described by a solid angle deficit in the (3 + l)-dimensional spacetime 
whose line element is given by 0] 

dv^ 

ds 2 = -dt 2 + + r 2 (d6 2 + sin 2 Odcf?) , (2) 

where the parameter a, which codify the presence of this object is smaller 
than unity. The energy-momentum tensor associated with this monopole has 
a diagonal form and reads: Tq = Tl = (a 2 — l)/r 2 and T 2 = T| = 0. 

Because we are interested in obtaining the thermal Green function, it is 
convenient to work in the Euclidean continuation of the Green function by 
performing a Wick rotation. So we shall work on the Euclidean version of 
the monopole metric above, expressed by the line element 

ds 2 = dr 2 + ^- + r 2 [dQ 2 + sin 2 Odtf) . (3) 
cr 

This paper is organized as follows. In section 2 we derive the Euclidean 
thermal Green function associated with a massless left-handed spin- 1/2 field 
in the pointlike global monopole spacetime adopting the imaginary-time ap- 
proach and using the Schwinger-De Witt formalism. In section 3 we 
present a formal expression for the thermal average of the zero-zero compo- 
nent of the energy-momentum tensor, {Too(x))p. Because this term cannot 



2 



be written in a closed form, its dependence on the temperature is not evi- 
dent. So, in order to obtain some quantitative information about its behavior 
we derive more workable expressions for specific values of the parameter a. 
Moreover, the complete information about its dependence on the tempera- 
ture requires that we proceed a numerical evaluation of some non-analytical 
contributions. Here in this paper, we develop these analysis in the high- 
temperature regime. In section 4 we summarize our results and give our con- 
clusion about the thermal bath of massless fermions in the global monopole 
spacetime. In the Appendix we derive explicitly the thermal average for the 
energy density for the special cases where a = 1/2 and give the dominant 
thermal contributions for the cases a = 1/3 and a = 1/4. 



2. Thermal Spinor Green Function 

The Green function associated with a massless left-handed spin- 1/2 field in 
the global monopole spacetime has been obtained recently M. There this 
Green function is presented by applying the respective Dirac operator on the 
two-components bispinor as we shall summarize below. 

In order to write down the Dirac operator in the spacetime described by 
(0) we choose the following basis tetrad: 



'(«) 



/ 1 \ 

a sin 9 cos (p cos 9 cos ip/r — sin ip/r sin 9 

a sin 9 sin ip cos 9 simp /r cos ip/r sin 9 

V a cos 9 — sin 9/r / 



(4) 



With this choice the Dirac operator for a massless left-handed fermionic field 
in a global monopole spacetime given by 



reduces itself toQ 



d t — aa r d r a d e 

r 



1 



r sin 9 



+ 



a 



-a 



(5) 



(6) 



1 Here = (I,a k ), a h , with k = 1,2,3, being the 2x2 Pauli matrices and T M the 
spin connection. 
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with o u = a ■ u, where u denotes the standard unit vector along the three 
spatial directions in spherical coordinates. 

The Feynman two-components propagator obeys the differential equation 

f L S F (x,x') = -L=8W(x-x')Iw (7) 

and can be given in terms of the bispinor G F (x,x') by 

S F (x,x') = if L G F (x,x') , (8) 
where G F (x,x') obeys the 2x2 differential equation 

CG F (x,x') = --±=5 {A \x - x')I {2) , (9) 

with 

C=u--R. (10) 
4 

In the above equation, the generalized d'Alembertian operator is expressed 
by 

□ = g^V,V v = gT (d„V v + T,V V - T« a V a ) , (11) 

and R, the scalar curvature, is equal to 2(1 — a 2 )/r 2 for this spacetime. 

From now on we shall work with the Euclidean version of the spinor Green 
function. The Euclidean bispinor, Ge(x,x'), can also be obtained from the 
Schwinger-DeWitt formalism as follows: 

G ( x,x) = J dsK{x,x';s) , (12) 

where the heat kernel, K(x,x';s), can be expressed in terms of the eigen- 
functions of the operator £ as 

K(x, x'; s)=Yl ^x(x)^ x (x')e- sX2 , (13) 

A 

A 2 being the corresponding positively defined eigenvalue. Writing 

C$( x ) = -A 2 $(x) , (14) 
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the complete set of normalized solutions of the above equation is given by 



A 2 



ap 



2nr 
E 2 + a 2 p 2 
/ + 1 1 



I 



a 



~2 ' V2 



1 

a + 2 



(15) 



(16) 



(k) 

where J u is the Bessel function of the first kind and f)^, with k = 1, 2, are 

the spinor spherical harmonics eigenfunctions of the operators L 2 and a ■ L. 
(The explicit expressions for these spherical harmonics and their eigenvalue 
equations are given in |4]]. There, it is explicitly shown the dependence of j 
with I.) 

]) we get: 



Now substituting (p~5|) into 



-"-OO ( X , X , Sj 



/OO fit.] roo 



sA 2 



1 



Aa\fnrr' s i / 2 



exp 



(AT) 2 q 2 +r 2 +r /2 



E 



2a 2 s 



_rr_ \ r (2) 



2a 2 s 



41 



(17) 



where J v is the modified Bessel function and Cj^. (fi, fi') = <£>}^. (O'). 
With the help of || we can see that for a — 1, V\ = v 2 = 1 + 1/2, it is possible 
to proceed the two summations on the quantum numbers j and rrij and (|T7|) 
reduces itself to: 



„(*) 



,(*)t 



x,x : s) 



'(2) 



(18) 



Now we are in position to obtain the Euclidean bispinor Ge{x,x') by 



2 In Ref. jl| we have used another parametrization for the time component. There 
goo = ot 2 ■ For this reason there is a small changing in the expressions obtained for the 
normalization of the field itself and its eigenvalue. 
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using (0). Again with the help of || we obtain: 

g e (x,x') = £ [g, 1 _ 1/2 ( M )c« j (^^') + 

Q„ 2 _ 1/2 («)cS.(^cr)] , (19) 

Q// being the Legendre function and 

a 2 (r — t') 2 + r 2 4- r' 2 

2rr' v y 

The Euclidean thermal bispinor can be obtained in the same way by using 
the thermal heat kernel, Kp(x,x'] s): 

POD 

Ga(x,x')= dsKg(x,x';s) . (21) 
Jo 

For an ultrastatic spacetime, i.e., static with g 00 = 1, Braden has 
proved that the thermal heat kernel can be expressed in terms of the sum 

oo 

Kf,{x, x'- s)= J2 (-l) n ^oo(x, x' - nXP; s) , (22) 

n=— oo 

where A is the "Euclidean" time unit vector. The zero-temperature heat 
kernel, K^ix^x'] s), given previously by (|17D can be factorized as: 

K^x.x'-.s) = K {1) (t,t';s)K {3) (x,x';s) , (23) 

where K^(t, t'\ s) and K@)(x, x'\ s) obey, respectively, the differential equa- 
tions below: 

(l-^) K ^^° < 24 > 

and 

" + i^) ^ (3)(f ' f ' ; S) = ° ' (25) 
Identifying in (fi7| ) as the solution for (El) 

Ar 2 

tf (1) ( r , r ';s) = l^, (26) 
V s 
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we can see that only K^(t, t'\ s) will be affected by the temperature, so we 
obtain 

Ko.)p(r, r'; s) = -= £ {-1)^—^ , (27) 



n=— oo 



which is antiperiodic in coordinate r. Now we can write down our result for 
the thermal heat kernel: 

1 1 00 
K p {x,x'-s) = - — i==^e E (-l) n exp" 



E 



(28) 



Finally substituting fl28|) into (|2T|) , we obtain an explicit expression for 
the thermal bispinor, Gplx, x'): 



1 oo 

+g, 2 -l/2K /3 )Cj 2 ^ i J (^^ , )] , (29) 



where 

a 2 (Ar + n(3) 2 + r 2 + r' 2 
u ^ = ^ ' ( } 

As we can see form (p9|), the thermal bispinor contains explicitly the zero- 
temperature Green function, Ge(x,x'), plus purely thermal contributions. 
Moreover, because the Legendre functions, Q^(z), decrease proportionally 
with for large arguments ||, we conclude that, for a given set of values 

of the quantum numbers j and rrij, the most relevant thermal corrections to 
Gp(x, x') above, come from the first terms in the series in n. 

Finally we can write down our Euclidean thermal spinor Green function 
substituting (|29D and the Euclidean version of @ into (§): 

Spix, x') = i (-id T - aa {r) d r + -a^a ■ L + G p (x, x') . (31) 
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3. Thermal Average of the Energy- Momentum 
Tensor 



The energy-momentum tensor associated with spin- 1/2 fields, T^(x), is bilin- 
ear function of the fields. Consequently its vacuum expectation value (VEV), 
(T^ u (x)), can be evaluated by the standard method using the Green function 

i- 

In M the renormalized VEV of the energy-momentum tensor associated 
with a left-handed spin-1/2 field in a global monopole spacetime was explic- 
itly calculated. There, it was shown that it has the following structure: 

CW>** = K + S XH] , (32) 

where the scaling parameter /i appears after the renormalization procedure. 
In fact, applying the point-splitting renormalization procedure and extract- 
ing from the spinor Green function all the divergences by subtracting the 
correspondent Hadamard function, the renormalized VEV of the energy- 
momentum tensor presents a scale-dependent logarithmic term. This is a 
consequence of an ambiguity in the definition of (T^)^ 9 ' in a even dimen- 
sional curved spacetime. The back-reaction probelm consists of solving the 
Einstein equaion in presence of quantum corrections. So, this one-loop equa- 
tion contains besides (T^)^ 9 ' as source of energy-momentum tensor, extra 
purely geometrical terms |7j . The coefficients of these terms must depend on 
the renormalization scale since the theory must be independent of this scale 
due to renormalization group equation. The tensors A 1 ^ and depend only 
on the metric parameter a. Moreover, because of the spherical symmetry 
these tensors are diagonal with components Ag = A 1 ^ and B d = Bf. There- 
fore there are only six unknown components. The renormalized VEV of the 
energy-momentum tensor must be conserved, i.e., 

TO R : 9 - = , (33) 

and gives the correct conformal anomaly M, which for massless spinor two- 
components field, reads || 

W>"* = i^™" = ' (34) 



8 



Here we have introduced the new definition T := rAt ™ 2 . The general form 
of the coefficient 02 may be found in [10 , which for the global monopole 
spacetime has the following form: 

1 — a 4 

a2 = _ ~60H~ /(2) ' (35) 



Taking into account ([33|) and (|3J) it is possible to express the tensors A 1 ^ 
and B v in terms of their zero-zero components and the trace T by 



A; = diag (a° , -T + A° + B° , T-A° 



B° B°\ 
^,T-A° -^) , (36) 



and 

^ = 5^(1,1,-1,-1) . (37) 

So the only problem left is to determine the components A® and B®. Using 
the point-splitting approach, the VEV of the energy-momentum tensor for a 
massless spinor field is given by: 

(T^(x)) = - lim tr [a^Vu - V u ,) + a u (V^ - V M 0] S F (x, x') , (38) 

by which we have 

(T 00 (x)) = % lim d 2 t tr (G(x, x')) = - lim d 2 T tr (G E (x, x')) . (39) 

x'— >X X 1 — >x 

It was pointed out in that only the time derivative in the bispinor gives 
a nonzero contribution to the zero- zero component of the energy-momentum 
tensor. Moreover taking the coincidence limit Q = Q' into (|T9"D, summing 
over rrij and using the integral representation for the Legendre function U, 
it is possible to develop the sum over j. Finally we arrive at the following 
expression for the Euclidean bispinor: 

I f°° dx 1 
^ (Ar ' r ' r ' )= 16^vA ^-^ sinh^^ J /(2); (40) 

where the function b is expressed in terms of the one-half of the square of 
the geodesic distance a in (r, r) surface 

b 2 = a 41 

2rr' y ' 
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with 

Ar 2 

2a = At 2 + . (42) 
or 

The Green function ((40|) is divergent in the coincidence limit b ^ 0. So, 
in order to obtain a finite and well defined result for (|39D , we must renor- 
malize it subtracting all its divergences. There, it was adopted the point- 
splitting renormalization procedure and subtracted from the Green function 
the Hadamard one, Gh(x,x'), which is expressed in terms of the square 
of geodesic distance a. The latter presents short distance behaviors pro- 
portional to 1/er, as a free theory, plus a logarithmic scale-dependent term 
proportional to the scalar curvature of the spacetime i?ln(/i 2 a). 

The thermal average of the energy-momentum tensor of a thermal bath 
of massless fermions in the global monopole spacetime can be computed 
in similar way replacing the ordinary spinor Greens function in ([38]) by its 
thermal correspondent, Sp(x,x'), given in fl3~l"|). For the sake of simplicity let 



us consider in this paper the zero- zero component of the thermal average of 
the energy-momentum tensor, which is also given by 

(T 00 (x)) p = - lim d 2 T tr (G p (x, x')) . (43) 

as'— >x 

The above expression gives us a divergent result which comes exclusively 
from the the zero-temperature contribution of the thermal bispinor, the com- 
ponent n = of (2~9), which we denominate by G^x, x'). In fact this singu- 



larity is consequence of the evaluation of the Legendre function at unity. So 
in order to obtain a finite and well defined result for (|43"D we have to renormal- 
ize its zero-temperature part only, (T o(x)) 00 . In this way the renormalized 
thermal average of the zero-zero component of the energy-momentum tensor 
can be expressed by 

(TM)f 9 = <T o(*)>r- + (Too(x))p . (44) 

Because the first term in the right hand side of the above equation has 
been computed before, in this paper we calculate the purely thermal correc- 
tion, (T Q0 (x))f3. Again taking the coincidence limit Q — Q' into (^J) and 
repeating the same procedure adopted for the zero-temperature case, we get 

(Tao(x))p = - lim d 2 T tr(G p (Ar,r)) , (45) 

Ar->0 
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where 



1 

167r 2 r 2 



E(-i) n r 



dx 



hi sinh 2 ( 



arcsinh(x) \ 

(2) > 



(46) 



being 6 2 = - ^rir^ ■ Because b n > 0, the purely thermal correction in 
( [44"D is finite. Changing the variable x by b n x in (flBf ), taking its second time 
derivative followed by the coincidence limit, we find: 



C?oo(aO) 



^2 r 4 



00 dzz 1 



n=1 h ^z 2 - 1 
za 2 nf3 cosh(y) 1 



2 cosh 2 (y) + 1 



sinh%) (l + 



Ar 2 



2r sinh (y) 



2 .2 Q ,2 n 2 / g2 
472 



3/2 



with 



arcsinh 



(47) 



(48) 



Unfortunately, from the above expression, it is not possible to obtain 
any information about the behavior of the energy density associated with 
a thermal bath of massless fermions in this spacetime as a function of the 
temperature. At this stage we only can observe that the integrand decreases 
exponentially for large values of z and that the series in n converge faster 
than 1/n 2 when n goes to infinity. Because we want to have a concrete 
information about the dependence of (Too)/? with the temperature, our next 
steps is to present this dependence for some specific values of the parameter 
a. In what follows we present this information for three different situations: 
a) a close to unity, small solid angle deficit (excess), b) small value of a, 
large solid angle deficit and c) large value of a, large solid angle excess. (In 
the Appendix we also present some concrete results for a = 1/2 and brief 
discussion for the cases a = 1/2 and a = 1/4.) 



3..1 Parameter a close to the unity. 

In this subsection we present the explicit dependence of the purely thermal 
correction to (To Q (x))/3 admitting that |1 — a>\ « 1. In fact, because the 
global monopole is a very heavy object which appears as consequence of a 
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breakdown global 0(3) gauge symmetry to U(l), for a typical grand unified 
theory 1 — a 2 ~ 10~ 5 . So in this case we can expand (47) in powers of the 



parameter rj = 1 — a . After long calculation we arrive, up to the first order 
in rj 2 , at: 



(Tao(x))p 



7tt 2 



120/3 4 

J2 oo 



1 + 2 V 2 ) + 



rj 



4vr 2 r 4 



E 

71=1 



-1^ 



arctanh ( 



n 



2r] 2 



7r 2 r 4 



E 

n=l 
A 



:-v, 



71° 



{l3/t) 



(f3/r) 2 ^l + n 2 f3 2 /Ar 2 
dz 1 

arcsinh(zn/3/2r) x 



6(f^) +19(f^) +12 



(1 + z 2 n 2 (5 2 /Ar 



2^3/2 



(49) 



Again this result also does not clarify completely our objective because 
of the presence of two contributions which cannot be expressed in terms any 
functions found in literature. However the above expression has some advan- 
tage when compared with the previous one obtained without approximation, 
Eqs. ( f4"T|) and (f48|). From (f49[) it is possible to proceed a numerical evalua- 
tion of both contributions for some specific values of the ratio £ := /3/r and 
in this way to provide a quantitative information for (T o)/3- Our numeri- 
cal analysis were developed in the high temperature or large distance limit, 
£ << 1, and the results are shown in Figs. 1(a) and 1(6) for «Si(£) and 82(C), 
representing, respectively, the first and second summation as shown below: 



5 (e) = y (~l) w 1 arCtanh Wl + n^/4 



n 2 e 



1 + n 2 £ 2 /4 



(50) 



and 



S»(0 



dz 



~ (-l) n 1 
6 (£|£) 4 + 19 (f«) 2 + 12 



r arcsinh(zn£/2) x 



(51) 



(l + ^ 2 n 2 £ 2 /4) 3/2 

From the graphs displayed in Fig. 2, which exhibits the logarithmic behavior 
of the previous summations, it is possible to infer the following dependence 
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for Si and S2 with £: 

= ^ (52) 

and 

&(0 = f 2 (53) 

where a x = -4.701 with q x = 2.228 ± 0.029 and a 2 = -8.735 with 
q2 = 4.000 ± 0.000. From these two results we can see that the most rel- 
evant contribution comes from 5*2, which, according to the precision of our 
numerical analysis, provides to ( |49"D a term proportional to 1//3 A , so indepen- 
dent of the distance from the point to the global monopole. 



3.. 2 Small Parameter: a << 1 

In the limit a << 1 the leading term in the expansion of the integrand of 



(|47D give us 



] 00 POO 



z 2 2cosh 2 (zn/3/2r) + 1 



V^^l sinh 4 (zn(3/2r) 



(54) 



which becomes independent on the parameter a. Once more this result does 
not elucidate our main question, which is to provide the information about 
the dependence of the thermal correction (T 0Q (x))j3 with the temperature. 
Again only numerical evaluation provides this information. The summation 
part, denominate by S3, is clearly a function of the ratio £ = (3/r. Again our 
numerical results in the high temperature regime, £ << 1, are exhibited in 
Fig. 3(a) for S3 itself and in Fig. 3(6) for its logarithmic. From both it is 
possible to conclude that 

St® = ^ , (55) 

where = —9.638 with 03 = 4.000 ± 0.000. This result indicates that 
for large solid angle deficit, the high-temperature contribution to the energy 
density is uniform and independent on the parameter a, i.e., (Too)/? ~ l//? 4 . 
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3.. 3 Large Parameter: a » 1 

Now in the limit a » 1, the leading term in the expansion of (|7]) reads: 

dz 1 



(Toq(x))p 



a 2 



2 



n=l 



3 y v tl 

^ n 2 A ^z T ^l[\n(aznf3/r)}' 1 



2n 2 r 2 {3 

« (_1)» p & 1 1 

^ n 2 A v^^Ifln^n^/r)] 3 ] ' 1 J 

The above result presents two terms whose their sum cannot be evaluated 
analytically. They appear as the first terms of an expansion of a well defined 
expression. Consequently they are valid only in the case where the argument 
of the logarithmic never becomes equal to the unity. This means that, from 
the numerical point of view, a should assume values bigger than some critical 
one, which depends on the interval domain of the variable £, in such way the 
divergence is never attained. In order to proceed a numerical evaluation of 
both terms, let us write down both terms in a convenient form, absorbing 
the factor r 2 /j3 2 = l/£ 2 . Defining them by S 4 and S 5 they are given as: 

00 (— l) n 1 f°° dz 1 
^[ n 2 4^ Ji y/z 2 - 1 [In {azn£)\ 

and 

00 (—i) n l r°° dz 1 



n=l 



rr 



£ 2 ii Vi 2 ^! [In (azn£)Y 



In this way there appears the factor 1/r 4 multiplying both terms in (|56|) . 
Adopting for ^ = (3/r values in the same interval used in the other previous 
numerical analysis our results for the two summations can be expressed as 
shown below: 

*4(fl = ^ (59) 

and 

S 5 (0 = ^ , (60) 

where we exhibit in Tables 1 and 2 the exponents and respective coefficients 
for different values of a for S4 and S5, respectively. Unfortunately we could 
not find any relationship between these parameters with a, although we can 
notice no significant variations on them, i.e., the exponents, g 4 and q$, do not 
change appreciably when we change a. The dominant dependence of (T o)/3 
with a is explicitly given by the power a 2 in ([56]). 
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4. Concluding Remarks 



In this paper we have considered the thermal average of the energy density 
associated with a thermal bath of massless spin- 1/2 fermions in the back- 
ground of a pointlike global monopole. The reason why we have taken into 
account massless field is because in this case it is possible to express the spinor 
Green function in a closed form in terms of special functions. Moreover, we 
decided to consider only one helicity state, the left-handed one, associated 
with the neutrino's degrees of freedom. In order to develop this analysis 
we have derived a general expression for the respective Euclidean thermal 
spinor Green function, Sp(x,x'). Using this Green function we obtained a 
formal expression for the the zero-zero component of the thermal average 
of the energy-momentum tensor, (Too)/?, which contains a zero temperature 
term plus purely thermal correction. Because this correction is expressed 
in terms of not solvable integrals, its dependence on the temperature could 
not be evaluated explicitly. So, to get some quantitative information, we 
decided to analyze this term for specific values of the parameter a. In this 
way three formal expressions were obtained for special cases related in last 
section. Unfortunately our results were expressed in terms of infinite series 
not analytically evaluated; however from them it is possible to develop a nu- 
merical analysis. We developed these analysis in a high-temperature regime. 
Our final results enables us to evaluate qualitatively the thermal corrections, 
(Too) p. For the cases a close and much smaller than unity, the corrections 
are dominated by terms of order T 4 , consequently independent of the dis- 
tance from the point to the monopole. For the case a » 1, there appears 
an explicit dependence of this correction with the distance. 

The obtained results may have some applications in the early cosmology, 
where the temperature of the Universe was really high. We can see that 
the thermal contributions modify significantly the zero-temperature energy 
density in that epoch. 

A few years ago, Linet |TT] has obtained the Euclidean spinor thermal 
Green function for a massless spin-1/2 fermionic field in the spacetime of a 
global cosmic string. There, it was possible to express the thermal average 
of the energy-momentum tensor in a closed form. Linet has found that the 
high-temperature dominant contribution for energy density is proportional 
to T 4 and that the thermal bath of massless fermions is not perturbed by 
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the presence of the straight cosmic string 3 . In the pointlike global monopole 
spacetime we cannot say that. For two specific cases treated here the results 
for thermal energy-densities in the high-temperature regime present explicit 
dependence on the parameter a. Only the case where a « 1 does not 
present this dependence. 

The contribution due to the temperature in the renormalized vacuum 
expectation value of the square of the massless scalar field and its respective 
energy density in the global monopole spacetime, were recently calculated 
in fl2| using similar procedure adopted in this paper. There appears a non- 
vanishing contribution to the energy density, coming from the geometry of 
the spacetime, to the dominant 1//3 4 term. 



A Thermal- Energy Density for Specific Val- 
ues of a 

Here in this appendix we present the purely thermal correction for the energy 
density for the cases where the parameter a assume very specific values. In 
this first part we take a = 1/2. For this case it is possible to develop the 
integral which appears in (fiB). It is given by 



64ir 2 r z 



1- * arctanh 1 



(A.61) 

with b n = ■ From the above expression we get, after some calculations, 

the following result: 

tm / ^ 7tt 2 1 1 1 , /4vrr\ 

+ 64|v COSh (if) C ° SSech2 




3 The high-temperature dominant contribution given in the Linet's paper, (Tqq(x)) — 
, is precisely twice bigger than the flat one found in our result in ( ^9[) . The discrepancy 
between both results is because we considered in this paper just one helicity state for the 
fermionic field. 
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E(-!) n 



arctanh! 



oo arctanh I , 1 = = 

3 (A.62) 



256jrV,^ ' (i + n'/P/ier 1 )"* 

In the above expression appears two contributions which cannot be eval- 
uated analytically. Again if we are inclined to know how these terms behave 
as a function of the variable £ we have to proceed a numerical investigation. 
The explicit expressions writing in terms of £ are given below: 



n=l 

and 



arctanh ( — -, — = ) 



oo 



arctanh ( 



i 



««)=E(-ir— (A.64) 
n=i (1 + n 2 £ 2 /lQ) 1 

From our numerical analysis, developed in the high-temperature regime, 
we obtained the following results: 

m) = ^ ( A - 65 ) 

and 

Sr(0 = ^ , (A.66) 

where a 6 = -0.552 with q 6 = 0.184 ± 0.019 and a 7 = -0.552 with q 7 = 
0.184±0.019. So we can conclude that these two terms present the same order 
of correction to (T o)/3 and are not relevant. So the dominant contribution 
comes from the term independent of the distance from the point to the global 
monopole, followed by a term of order l/r 2 (3 2 . 

For the case a = 1/3 the procedure is similar to this previous one. The 



integration in (^) can be done and after some intermediate steps it is possible 
to conclude that the dominant contribution for the energy density, in high- 
temperature limit is also given by 



(Too(x))p = — -57 + hns 2 • (A.67) 
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The same dominant behavior is obtained for the case a = 1/4. We shall not 
repeat the intermediate steps. 

As our final comment in this appendix, we want to say that the domi- 
nant contributions for (Too(x))p in these three examples, coincides with the 
result found for a flat spacetime, (T 00 (x))/3 = . However, sub-dominant 
contributions are consequence of the non flat geometry of this spacetime. 
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B Figure Captions 



Figure 1: The figures 1(a) and 1(b) exhibit, respectively, the behavior of S± 
and S 2 with £ in the region [0.01; 0.1]. 

Figure 2: This figure exhibits the logarithmic behavior for both functions, 
Si and S 2 with the logarithm of £, in the region [0.01; 0.1]. From it, it is 
possible to estimate the leading dependence for both quantities with £. 

Figure 3: The figure 3(a) exhibits the behavior of S3 with £ in the region 
[0.01; 0.1]. The figure 3(b) exhibits the logarithmic behavior for S3 with the 
logarithm of £, in the region [0.01; 0.1]. From it, it is possible to estimate its 
leading dependence with £. 

Table 1: The table 1, exhibits the exponents g 4 with their respective coeffi- 
cients a 4 of £4 for different values of a. 

Table 2: The table 2, exhibits the exponents q§ with their respective coeffi- 
cients a 5 of S 5 for different values of a. 
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Table 1 


a 






1000 


2.930375608 ± 0.1669581160 


-2.596702947 


1500 


2.8414692776 ± 0.1335029194 


-2.348446259 


2000 


2.788466901 ±0.1156906043 


-2.193442147 


2500 


2.751914554 ± 0.1042923505 


-2.082895068 


3000 


2.724562303 ± 0.09622178978 


-1.997972468 


Table 1: 


Table 2 


a 


<?5 


a 5 


1000 


2.632887909 ± 0.12000918880 


-3.059963397 


1500 


2.569103696 ± 0.09581311868 


-2.887785837 


2000 


2.531144121 ± 0.08291154623 


-2.780842812 


2500 


2.505004736 ± 0.07464879774 


-2.704862735 


3000 


2.485469694 ± 0.06879530831 


-2.646669348 



Table 2: 
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